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PRECED,Nq p ** bunk not nu® 

ABSTRACT 

The reflection of short laser pulses froa the ocean surface is 
analyzed based on the specular point theory of scattering. The expressions 
for the averaged received signal, shot noise and speckle induced noise are 
derived for a direct detection system. It Is found that the reflected 
laser pulses have an average shape closely related to the probability 
density function associated with the surface profile. This result Is 
applied to estimate the aean sea level and Significant Wave Height from 
the receiver output of the laser altimeter. 
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1. INTRODUCTION 

In a previous paper by Gardnar [I] • cha statistical charactwrlatlca 
of short laaar pulaaa (pulaa length laaa than 1 cm) that have baan reflected 
from tha ground vara atudlad. In this report, wa axtand tha pravious 
rasults to tha caaa vhara tha raflactlona occur at tha ocaan surfaca instaad 
of froa a ground targat. 

At naar vortical incidence, tha raflactlon of laaar radiation froa tha 
ocaan is aalnly accounted for by tha scattering froa specular points which 
are randomly distributed on tha ocaan surface. If tha laser pulse length 
is short conpared to the surfaca height variations* tha raflactad pulses 
from the oc.'an will have an average shape related to tha height probability 
density of tha specular points as wall as tha overall probability density 
of tha height variations. Therefore: a short-pulse laaar altimeter can be 
used in tha determination of sea states aa wall as the mean sea level. 

Short-pulse satellite altimeters in the microwave range have been used 
with good rasults in observing sea states [2], [3]. For radar altimeters, 
the transmitted pulse width is typically a few nanoseconds or longer, end 
the antenna beam width is usually on the order of degrees. For laser 
altimeters, the transmitted pulse width can be as short as tens of pico- 
seconds, and the laser divergence angle can be as small as 10 urad. Because 
of these significant differences in the transmitter parameters, previous 
results oe the radar altimeter have to be examined before they are applied 
to the laser altimeter. 

In this report, we first derive the received signal in Chapter 2. 

In Chapter 3, ocean surface statistics are discussed. The temporal 
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moments of returned pulses srs examined In Chsptsr 4. In Chipter 5, 
wavs forms of ths received signal ere calculated uaing both Gauaslan and 
non-Gauss lan ocean surface statistics. It is shown that, with comparable 
parameters, our results for laser altimeters parallel those of radar 
altlmetera. Finally, in Chapter 6, the effect of non-normal incidence is 
considered. 
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2. RECEIVED SIGNAL 

The geometry of the laoar altimeter and ocaan aurfaca la llluatracad 
In Flgura 1. Tha altltuda of cha alclnaear measured from cha naan aaa 
laval la a. Tha 2-D aurfaca prof 11a la daacrlbad by $(£>, and lta 
corraaponding alopa In tha x- and y-dlractlona la danotad by £ (£> and 
respectively, vhara £ la tha horizontal position vac tor on tha 
ocaan aurfaca and la aaaaurad from tha cantar of tha laaar footprint. 
Initially, va aaauaa that tha laaar pulaa la incldant normal to tha maan 
aurfaca laval. In Chapter 6, the affacta of non-normal Incidence are 
dlacuaaad. 

Ua first derive tha mutual coherence function at tha receiver plana. 
Tha cooplex amplitude of the pulaad laser beam at the tranamlttlng 
telescope is 

lu> t 

u T (r,t) - f(t) a T (r) e 0 (1) 

where 


f(t) • laaar pulaa amplitude 

a^(r) ■ complex amplitude croaa-aaction of tha laaar beam 
u) q ■ laaar frequency 

r_ ■ (x,y) - transverse coordinate vector. 


Using tha Fresnel diffraction formula, tha field incident on tha ocean 
surface la 


B t (r,.,t) - T 2/2 f 


t “ 


2z 


2cz 


a 1 (r,z) axp 


. o 2 

ut-kz- — r 
o o 2z 


( 2 ) 
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Figure 1. Geometry of the laser altimeter and ocean surface for normal 
incidence. 



T la tha intensity transmit tr.uca of tha atmosphere, k la cha wav •number 
A o 

and X Q la tha wavelength of tha laaar radiation. Throughout thla papar, 
tha apatlal lntagrala ara assumed to be avaluatad over tha antlra plana. 

In darlvlng Equation (2), va have aaauaad that tha raa laaar pulaa width 
(Oj) and tha araa of tha transacting talaacopa (A^) aatlafy tha condition 

K 

co f » — (4) 


For naar normal incidence, tha raflactlon of laaar light from tha ocarn 

aurfaca la mainly accountad for by acattarlng from randomly dlatrlbutad 

apacular polnta on tha aurfaca. Slnca tha ocaan aurfaca la vary rough on tha 

acala of tha optical wavelength, tha acattarad powar will ba proportional to 

tha nuabar of apacular points lllumlnatad and tha acattarlng crosa-sectlons 

of apacular polnta. According to tha raaulta of Kodla [4], aach lllumlnatad 

apacular point acattara Ilka tha tangant aphare whoaa radlua la tha geometric 

maan of tha two principal radii of tha aurfaca at the apacular point. In 

the optical region, tha acattarlng cross section of a sphere with radlua R 
2 

Js Just ttR [5], Than, tha reflected field in tha plana Immediately above 
the ocean surface la given by 


U g (r,*,t) - R(0) [imjr^rj ] 


1/2 


2C(r) 


IT. 


r,z,t + 


(5) 


where 


R(0) ■ Fresnel reflection coefficient for normal incidence 
5(0 - surface elevation at jr 

| r^r^ | ■ absolute value of product of principal radii of curvature 
at the specular point 

n • number density of specular points at the point of reflection. 
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|r r. I can ba axpraaaad In tar na of tha darlvatlvaa of tha aurfaca prof 11a 
a d 1 

C(r) aa follow* [6] 



(1 * ♦ ^) 2 
Sy " 


1 5 , 


( 6 ) 


whara 






Equation (5) lapllaa that whan tha lncldant flald la raflactad It la 
raducad In anplltuda and undargoaa a propagation dalay which la proportional 
to tha aurfaca halght dlatrlbutlon. Tharafora, tha flald In tha plana of 
tha racalvlng talaacopa la 


lw Q t T a 

U(r,*,t) - a(r,z,t) a - — axp 

o 


o» t - 2k » - ^ r^ 
o o 2z 


| <* 2 £ 


R(0) (tm|r^r b | f 


c cc c 


>axp 


-ik 


£.*£ 




(7) 


In darlvlng Equation (7), tha laaar pulaa width and racaivar apartura 
araa ara aaaunad to aatlafy tha condition 


co f ” T 


( 8 ) 
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Th- itutual coherent’ function In defined ee 


< 9 ) 


where the angle bracket denotee tvpectatlon with reepect to epeckle end 
the nicroe true cure of the surface. 

Under th. assumption that the ocean surface is rough on the scale of 
the optical wavelength and chat the nicros true cure is unreaolvable by the 
receiving telescope, we can first perform the expectation over speckle, 
and Equation (9) becomes 


J a < -l ,t l ; -2 ,C 2* “ T a*~ 2 I *K°> I 2w «*P 


- 1 2! <*i 2 “ r 2> 


<n(£,0 


(i + il + e 2 ) 2 


^xx ^yy " ^xy 


2 > f - i^) f (t 2 - axp 


j d 2 £.|e l (fi.,e) 


^£(£i-r 2 ) 


( 10 ) 


who re 


* 



2 £<£) 


(ID 


The expectation inside the integral in Equation (10) Is with respect to 
the number of scacterers and their scattering cross- sect ions. By using the 
results of Barrlck [6] and taking into account the £ dependence, we have 


<n(£,0 


(1 


♦ £ 2 ♦ £ 2 ) 2 
J£ JL_ 


yy 


- 5 


*u ♦ c 2 ♦ c 2 ) 2 p(c x ,c y |c(£)) 


( 12 ) 


*y 


where p(£ ,£ |0 is the Joint probability density function of £ and £ 
x jr x y 

at a given elevation £. 

Since che ocean surface can be 'c-icribed by the equation 


f (x,y , z) - * - £(x,y) - 0 


(13) 
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the surface normal sc chs specular point la oV.alned by caking the gradient 
of the left side of Equation (13) 


n • Vf ■ i • i( x - yC y (1*) 

AAA ^ 

where x, y, z ere unit vectors In the x-, y- and z-dlrectlons, and n la the 
normal vector at the specular point. If we dec te the angle between the 
t-axls end local surface no real n as 8, then froe Eq. (14) end Figure 2 we 
find that 9 Is related eo the surface slopes at the specular point by 

tau 8 - ♦ ( 2 ) 1/2 (15) 


For bsckscatter, 9 le also the Incidence angle for the specular point. 
Equation (15) laplles that only those specular points which ere tilted 
so as to be noraal to the incident wave ere effective In contributing to 
the received field. 'Ails is not an unexpected result. Since 8 Is also the 
Incidence sngle, it can be related to the horizontal position of the specular 
point p by 


tan 8 ■ ^ (16) 

Note that we ere assuming that £ is neasurad froe the center of the laser 
footprint. By substituting Equations (12). (15) and (16) into Equation (10), 
we obtain 


■ ■ 

I a ( -l* t l ; -2 ,t 2 ) * T a*” 2 2 ! (r i “ r ]} |*<0) | 2 w j d 2 £ 


•I *> I 

r ^ - • 

••*P ft £*(l x - £ 2 ) 


p( Cx*CyUV ) f< W ^2‘P 5 


(17) 


Ue assuee that the receiver f ield-of-vlew end interference filter are 
adjusted so that all the signal energy which is incident on the telescope 
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objective Is focussed onto the photodetector. Therefore, the total 
racaivad signal powar Is glvan by 

P(t) - | d 2 r W(r) |a(r,a,t) | 2 (18) 

where W(r) la sn appropriata apartura weighting function. W(j r) la equal to 
one for i: Inside the aperture and zero otherwise. For direct detection, 
the mean and covariance of the signal at the receiver output can be calculated 
using Caatpbell's theorem [7] 


E [s(t>] - E[P(t)] * h(t) 

o 

a* 

C (t., tj) - jjJ- [ dx E[P(x) ] h(t. - t) h(t 2 - x) 

J " o * 1 

-a 


+ 



dx 2 c p ( T 1 * T 2^ h ( t ^ 


(19) 


t 2 ) 


( 20 ) 


r\ la the efficiency of the receiver optica and detector, hf Q la the energy 
of one algnal photon, h(t) la the impulse response of the receiver electronics, 
and 0^ la the temporal covariance of the received signal power. 

The variance of a(t) may be regarded as the signal induced noise. The 
first term In Equation (20) is the quantum or shot noise component of the 
covariance. The second term arises from the randomness of the scattering 
cross-section, surface profile, and the speckle Induced fluctuation. 

Since the amplitude of the received field for a given realization of 
the surface profile is a circular complex Causslan process, the expected 
value and covariance of P(t) can be written In terms of the mutual coherence 
function of the received optical field. 

E[P(t) | 5] • | d 2 r J a (r,t;r,t) W(r) (21) 



Ths unconditioned mean received power is obtelned by taking the expectation 
with respect to surface profile. 


E[P(t)] « E[E[P(t) | £] ] 

- A R T^s" 2 |R(0) 1 2 ir j d^Isj^Cfi.,*) I 2 
•|f(t - *)| 2 


2 

l+*r 


| d* p(£ x ,£ y ,£) 


( 22 ) 


where p(£ ,£ , £) is the joint probability density of £ . £ and £. 
x y x y 

By substituting Equation (22) into Equation (19), the taean wave fora 
at the receiver output can be expressed as 


E [ s ( C) ] - A^T 2 x” 2 1 R(0) | 2 it | d 2 £|a t (£,z) 


^ + 2 
z 


2z p 2 2^(£) 

« pu x .yo *(t -r - « + — > 


(23) 


where g(t) - |f(t)| * h(t) is the point target response of the laser 

altimeter. In Equation (23). the Integration over £ is recognized to be 
the convolution of altimeter point target response g with the joint 
probability density function p(£ , £ ,£). 

•• y 

The power covariance function of the received field is defined as 


C p (t l’ t 2 ) " E{P(t l ) P(t 2 )} ' E{P(t i )} E(P(t 2 ) ) (24) 

The first term in Equation (24) can be written in terms of the amplitude 
of the received field as 


E{P(t 1 ) P(t 2 )} 
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where Che inner engle bracket la the expectation with respect to the speckle 
and micros tructure of the surface, and the outside expectation Is with 
respect to the surface profile £• 

Using the properties of the circular complex Gaussian fields [8], the 
expectation over speckle and microstructure Is carried out and expressed 
in terms of mutual coherence function 

UKt^ P(e 2 )> - J d 2 ^ | d 2 r 2 W(r 2 ) w(r 2 ) St-V-rV^l’ 'l' • , a (l 2 ' t 2 i -2’ t 2 ) 

+ |J t (r l t 1 ;t 2 .t 2 )| 2 ] (26) 


The expectation inside the Integral is over the surface profile, which can 
be written explicitly in terms of the probability distribution functions 
of the surface profile. The first term in Equation (26) becomes 

| d 2 r x d 2 jr 2 WC^) W(r 2 ) E [ J a (E 1 ,t 1 ;£ 1 , tj) J A ^L 2 ,t 2 »-2 ,t 2^ 


- T 4 z" 4 |r(0)| 4 it 2 


4 f d2 % ( <, %l a i ( %- z) l 2 l |a i ( £ 2 - 2) 


2 ' 

°1 

1+ 4 

z 


2)2 

°2 

1 + -i 

z 




d5 2 p(5 x 1 ,5 y 1 ^l ) P ^x 2 ’S 2 ^ 2) p(5 1 ( ^1 ),5 2 (£ -2 )) 


If (t L - *j) I 2 |f (t 2 - * 2 ) | 2 


(27) 


where p (C^ (p 2 ) • is the J oint Probability density of surface heights 

at and £ 2 » 

The second term in Equation (26) can be simplified by noting that the 

mutual intensity function only depends on the difference coordinate (r^ - £ 2 ) . 

-1 + -2 

By making the change of variables £ ■ r ^ - jc 2 and r^ - In Equation (26), 

the integration over r can be carried out to give 
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{ { *% "(ij) «(I 2 ) «|J.(» l .t 1 »£j.t 2 )| 2 l 

- I i 2 r ^(r) *[|.J l (r/2.t 1 i-r/2,t 2 )| 2 ] (28) 

where Ry la cha autocorralaelon function of tha racalvar apartura 

Ry (r) - | d 2 £ W(£) W(£ + r) (29) 

Notica that tha maximum valua of Ry occura whan .r • 0 and la equal 
to the aperture area. Normally, tha diameter of tha receiver apartura will 
be large compared to tha apackle correlation length [9]. In thia case, 
Ry(r) will be approximately conatant over tha important area of integration 
so that Equation (28) can be approximated by 

d 2 r Ry(0) E[J a (r/2,t i; -r/2,t 2 )| 2 ] 

- R w (0) | d 2 r E[|J a (r/2,t 1 ;-r/2,t 2 )| 2 ] (30) 

The integration over r in Equation (30) can be eaaily evaluated. By taking 
the expectation over £, we have 

R y (0) | d 2 r E[|J a (r/2,t i; -r/2,t 2 )| 2 ] 

2)4 r 

- X 2 tV 2 |r( 0 )| 4 d 2 £ |a (g.,*) I 4 1 + ^ dC P 2 (C X ,C Jo P(C) 

* 7 

•|f(t 1 - *)| 2 |f(t 2 - *)| 2 (31) 

Finally, the power covariance function of the received field can be 
obtained by aubatltutlng Equationa (22), (27) and (31) into Equation (24) 
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-4 

t 


yw 1 x o T i «' J i R «» i 4 *\ | <i £ i* t (£.*)r 

•j « P 2 (5 x .« j ,I« P<£> Ifttj - *)| 2 |f(t 2 - *)| 2 
+ I*«- 4 |R(0)| 4 p 2 * 2 | d 2 £l | d £ 2 l« 1 (£ 1 .P)| 2 |t 1 (£ 2 , I )| 2 

I d? 2 | dc 2 |f( tl - p^l 2 |f(t 2 - * 2 )| 


2 

l+! I 

z 


l + i 

s 


* [pte* .5 y U,) P<£ .£„ |c 2 ) p(£,.£ 2 ) - p(£ .£ ,£.) p(£ »£ ,£,)] 

*1 A 2 y 2 4 A * *1 1 *2 *2 C 

(32) 

The covariance of the output signal la obtained by substituting Equations (22) 
and (32) into Equation (20) 

2 


C s <t l- t 2 > " 


hf 


TV 2 |R(0)| 2 irA* f d 2 £|a 1 (fi.,£)| 2 


2 

1+£ 2 

z 


OD 

1 1 pU x .lyO | dt|f (t - *)| 2 h(t L - T) h(t 2 - r) 


hf 0 


i2 


JtV 2 |R(0)| 4 n\ f d 2 £|a t (£,«)!* 


2 

z 


dC 


*P 2 (5 x ,C JO p(0 g(t x - *) g(t 2 - m) 


t2 


hf 


T 4 z 4 1 R(0) | 4 w 2 ^ j d 2 £ x | d 2 p 2 |a 1 (£ 1 ,z)| 2 |a 1 (£ 2 ,z)| 


\ 

2 

2 

f 2l 

2 


1+! I 


°2 

1 

f M 

z 

. j 


z 

k 4 



d *l { d *2 8(t l “ V 8 ( t 2 " 


* lp( 5 x 1 ,C y 1 ^l ) p( 5 x 2 »S 2 ^ 2 ) p( W " p(C x 1 ’S 1 ’ 5 l ) 


P ( • £ v » O) 1 
*2 y 2 * 


( 33 ) 


15 


3. OCEAN SURFACE STATISTICS 


The wave height at Any given point on the oceen surface la the resultent 
of many wave components that have been generated by the wind in different 
regions and have propagated to the point of observation. Since the inter- 
actions between each wave component are weak [10], their notions are aasuaed 
to be weakly correlated. Therefore, under the central- Unit theorem, we 
expect the distribution of wave height to approach a Gaussian. 


The first approximation to the distribution of wave height is then 


p(C) - Uito 2 )” 1 ^ 2 exp 


_ JL 


(34) 


where o^ is the ras wave height. 

Significant Wave Height, SVH, is defined as the average of heights 
(fron crest to trough) of the one-third-hlghest waves observed at a point. 
It Is approximately equal to 4 times the rms wave height of the ocean 
surface [3]. 


SWH - 4 (35) 

The rms wave height can be related to the wind speed by Integrating 
the Phillips’ wind-wave-height spectrum [11] 

■ 0.016 W 2 meter (36) 

where W Is the averaged wind speed In meters per second measured at 
12.5 meters above the sea level. 

The first approximation to the joint distribution of the wave slopes 


is also Gaussian [12] 
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P(5 X .C y ) 



2 2 
i + 5 
x 


(37) 


2 

where S la the naan square valua of cha total alopa, daflnad aa 

S 2 - <e 2 > + <5 y > (38) 

2 

An amplrlcal ralatlonahip batvaan S and wind apaad la glvan by Cox and 
Hunk [13] aa 

S 2 - 0.003 + 0.00512 W (39) 


where W la in meters par aacond and S la dimensionless. However, the actual 
surface profile cannot be an exact Gauaalan for two reaaona. First, the 
wave height can never go to Infinity. Second, due to the presence of weak 
nonlinear Interactions between wave components, the actual ocean surface la 
not symmetric about the mean, as predicted by a Gauaalan distribution. This 
la born out by the observation that, on the ocean surface, wave crests tend 
to be relatively high and sharp, while the wave troughs are comparatively 
smooth and shallow. Mathematically, this fact could be taken Into account by 
including higher-order terms containing skewness and kurtosls In the probability 
density functions given by Equations (34) and (37). 

It is simple to show that surface elevation and surface slopes are 
uncorrelated at the same ?olnt. Therefore, If a Gaussian surface profile 
is assumed, the surface elevation and slope are independent. Thus, the 
joint density function of surface height and slopes can be factored as 


p(e,e x .e y ) • p(e x .e y > p(0 


( 40 ) 
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One important Implication of Equation (40) la that tha acat taring croaa- 
aactlon will ba tha aama from wava craat to vava trough, Indapandant of C. 

But, tha experimental raaulta obtalnad by Yaplaa at al. [14] for a 1 CHa, 

1-na pula ad radar Indlcatad an approximately llnaar Incraaaa of acattarlng 
croaa-aactlon from wava craat to vava trough. On tha othar hand, axparlmantal 
raaulta obtalnad In tha optical ranga, ualng N 2 » YAG and C0 2 laaara, lndlcata 
that maximum raf lactloa may occur at tha vava craat [IS]. Tharafora, tha maan 
aaa laval "aaen" by tha altlmatar would ba dlffarant from tha trua aaa laval. 

Jackaon [16] darlvad tha axpraaalon for tha joint probability danaity 
function of vava halght and alopa, which takaa Into account tha non-Cauaalan 
bahavlor of tha aaa vavaa. Hla raault la tha following 


pU x .O 


[2.(o* 


exp 


2 n 


i! + Js_ 

°l <> 


X 3 

l + r 


Q 

1 

■ 

3 

0 

V 

a. 

_ 2 





c 

1 

A 

K 

u/ 

V 



(41) 


whara la tha akavnaaa coefficient which la dafined by 



In arriving at Equation (41) , Jackaon aasumad a long craatad or 


corrugatad aaa aurfaca and mada uaa of tha one-dlmanaional Philllpa' 
aaturatad wava number apactrum. However, Equation (41) doaa not include 
tha affect of capillary vavaa, vhoae atatlatlca are a till uncertain. 
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Capillar, vavaa art wavaa with vavalangtha leaa than approximately 1.7 cm. 
Thay will not ba aaan by a radar altimeter with wavalengtha of 10 ca or 
larger, but will ba aaan by tha laaar altlaatar. Therefore, tha validity of 
Equation (41) In tha optical ranga raaalna quaatlonabla. In tha abaanca of 
a aora accurata axpraaalon, wa aaka uaa of Equation (39) In aoaa derlvatlone, 
ka aping In Bind that tha actual akawnaaa taraa could ba dlffarant In both 
aagnltuda and algn. 

For vartlcal Incidence, va naad tha axpraaalon for p(£ ,( ,£) avaluatad 

^ y 

at £ ■ £ • 0. Although tha axpraaalon for p(£ ,£) alona la not a no ugh for 

^ J ^ 

ua to obtain tha ganaral axpraaalon for p(£ ,£ , £), It doaa allow ua to 

^ y 

2 2 

obtain tha daalrad axpraaalon for p(0,0,£). By aubatltuting S for <£^> 

In Equation (41) t and evaluating at £ x - 0, v« have 


p(0,0,£) 


S 2 o*]“* 1 axp 


li! 
2 2 


3 

1 + r 


£ 

2 

°C 



(43) 


Equation (43) la vary almllar to tha marginal dlatrlbutlon of wave height 
with tha akawnaaa term included 


p(£) 







(44) 


It will ba ahown that. If Equation (43) la uaed In Equation (23) for tha 
returned algnal, the mean aea level aaan by tha altimeter la blaaed by 
tha amount A^o^. With the Inherent high accuracy of laaar a^timetera, thla 
biaa la elgnlf leant, and ahould ba conaldarad In tha proceaelng of received 

algnala. 
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4. TEMPORAL MOMENTS OF THE RECEIVED SICNAL 


In chin chapter, temporal aoiuei of the received signal art darlvad. 

For simplicity, w<t assume tha surfaca profila la Gaussian, so chat we 
can usa Equation (40) to factor tha joint dansity of surfaca slopes and 
alavacion. 

Tha temporal aoaencs of tha racalvad signal ara ralatad to tha statistics 
of tha ocaan surfaca profila. Tha k th -order aoaant is daflnad aa 

dt t k S(t) (45) 

Tha zeroth-order aoaant, m Q , is proportional to tha total racalvad 
signal anargy. Tha expected valua of a Q is calculated using Equations (23) 
and (40) 
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<N> is tha axpactad nuaber of datactad signal photons par racalvad pulsa, 

G Is tha gain of tha receiver electronics, Q is tha total anargy trensaltted 
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p«r lutr pula*, and la th* aqulvalant powar raflactlon coafflclant of 
cha anclra laaar footprint. Ganarally, B f la a function of tha aaa atata 
within tha laaar footprint. 

To coaputa tha atatlatlca of tha taaporal aoaanta, It la convaniant 
to axpraaa cha naan and covarlanca of tha racalvar output In tana of 
<« o >. Proa Equation* (23), (33) and (46), wa obtain 

E(a(t)] - <K> | d 2 £.b 2 (£,a) J d* PU) g(t - *) (51) 
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For n ■ 2, b 2 (£,t) la cha normalised affective croas-secclon of cha laaar 
beam. K ( la cha raclo of cha receiver aparcura araa co cha apackla 
correlation araa. Tha apackla correlation araa la Inversely propor clonal 
co cha araa of cha laaar foocprlnc. Typically, la such graacar chan 
ona [9]. 

Va ara primarily in tar •• cad in cha nornallsad moments 
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dc t K S(C)/ 
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If va assume cha fluccuaclona in S(c) dua co shoe noiaa and apackla 
ara small, cha maan of m^/m o can be written as [1] 
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Tha variance of a is 
o 

Var(m o ) - j" f dc 2 C # (t, .tj) 

— m —m 

By subsdCudng Equation (52) into Equation (58) and carrying out tha 
int«gradona , we obtain 


(58) 
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o s 


(59) 


Tho first-order normalised aomenc is cha time delay baevaan tha firing 
of tha laaar pulse and cha cancroid of cha received pulse. This delay 
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can b« uaad co eatlnata cha altitude of cha altimeter above cha taa level. 
If ve denote the received signal cine delay by T # , then the expected delay 
can be written aa 


<T > • E 

e 



(60) 


By using Equation (54), Equation (56) can be evaluated to give 

<V - j * 2 £ b 2 te.*) J d* P(0 p (61) 

Using the expression in Equation (9) for p, Equation (59) can be slw>Ufi«d 
_2 

vr\ 

(62) 


<T > .1** + 
s c cz c 


where 


\ m | d 2 £ X) 
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For Cauaslan surface profile, <£> ■ 0. It will be shown later chat 
if the non-Gauss lan nature of the surface profile is taken into account, a 
tern proportional to the skewness coefficient Is added to the right side of 
Equation (62). The expected delay given In Equation (62) la composed of two 
terns, z is the altitude of cha altimeter neasured from neon sea level. 


Therefore, — represents the round- trip propagation delay along a line normal 
c 

to the surface. la the rms wld£h of the effective laser footprint seen by 

the receiver aperture. The tern — represents the additional propagation delay 

cz 

resulting from the curvature of the diverging laser bean. 

The mean-square width of the received pulse la 


j 2 • J" dt(t - T g ) 2 S(t)/f dt S(t) 


( 64 ) 


The expected value of the pulse width can be calculated to give 
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E[o*} " + °f + "2 + (C,) ’ 2 d2£ - (p2 " °r )2 b 2^* r) (65) 

c 

2 2 
where cr^ ie the mean-square width of h(t), a f is the mean-square width of 

2 2 

the transmitted letter pulse ( | f | ) , end o* Is the verlence of the oceen 

2 2 

surfece profile, o^ end in Equetion (65) cen be computed from the known 
peremeters of the altimeter. Both the third term end the lest term depend 
on the see states within the footprint. Therefore, the width of the received 
pulse cen be used to obtain Information about ocean surfece conditions. 

This is discussed in detail in the next chapter. 
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5. WAVEFORM OF THE RECEIVED SIGNAL 


The later cross-section and waveform of the tranamlttad laaar pulaa 
are at turned to be Gauaalan in ahapa 




axp 


2o 2 (t) 
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where E q la the total energy of the transmitted later pulse, and o^(z) is 
the rms width of the transmitting laser cross section. o^(z) Is related 
to the divergence angle of the transmitted laser beam 9^ and the altitude 
of the altimeter z by 


o. ■ z can 0_ 
1 T 


( 68 ) 


In addition, the response of the receiver is assumed to be given by 


h(t) 


G(2wo 2 )' 1/2 




(69) 


where G is the gain of the receiver. 

The point target or Impulse response of the altimeter, defined 

2 

previously as g(t) ■ |f(t)| * h(t), can be expressed explicitly by 

using Equations (65) and (67) 
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If Gauaalan aurface aeaelatlca art aaaumed, we can aubatltuca 
Equation! (34), (37), (66) and (70) into Equation (23) and carry out 
tha integration. The raault la 
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1 - arf 
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and erf(*) ia tha error function. 


In arriving at Equation (72), we have made uae of the fact that 


2 ' 

1+£ 2 

z 


- 1 within tha laaer footprint. Tha juatif lcation behind this 

v " p 2 2 

la that the maximum value of la tan 9_, and 9_ la on the order of 

2 T T 

-3 * 

10 radiana or smaller for a typical laaer altimeter. 

Equation (72) ia found to be of the aaae form aa the reaulta obtained 

by Fedor et al. [3] and Hammond at al. [17] for radar altlmetera. 

Equation (72) la plotted in Fig. 3 through Fig. 7 for different acts 

of altimeter parametera and varloua aea atatea. It la obaerved that, for 

a beam divergence angle 10 -2 radiana or larger, the received waveform la 



(1) 10* rod 

(2) 3X10* rod 
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Figure 3. Mean received waveform* for laser altimeters with 
different beam divergence angles. 







RELATIVE AMPLITUDE 



Figure 5. Mean received waveforms of Che laser aldneter for 
different wind speed and SWH. 
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Figure 6. Mean received waveforms of Che laser altimeter for 
different wind speed and SWK. 
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Figure 7. Mean received waveforms of the laser altimeter 
for different wind speed and SWH. 
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highly asymmetrical, and tha figures obtained are similar to those obtained 
for radar altimeters [3], [17]. In this ease, the rise time of the leading 
edge of the received pulse can be used to Infer the roughness of the ocean 
surface, but the trailing edga of the received pulse Is relatively lnaen- 
sltlve to sea states [see Ref. 3]. 

For a beam divergence angle 10 radians or smaller, which Is typical 
for a laser altimeter, the received waveform Is found to be nearly 
Gaussian In shape. In this case, the centroid and rms width of the 
received pulse are easily Identified. The expected delay, given In 
Equation (62), Is 

<T > . 2* + ^tan“ 2 e T + 2 S * 2 )" 1 ( 76 ) 

see T 

where S le the mean square value of tha total slopes, defined previously 

In Equation (37) . Equation (76) differs from the previous result for 

_2 

reflections from the ground target [1] by the presence of the S term. 

Ground reflections ere primarily diffuse, end the reflectivity of the 

target is more or less uniform within the laser footprint. On the other 

hand, reflections from the water surface depend on the occurrence of 

specular points which satisfy the required slopes. From Equations (15) end 

(16) , we see that the surface slopes required for the contributing specular 

points ere larger for points further away from the center of the laser 

footprint. Since the occurrence of specular points with large surface 

slopes Is less probable than that with small surface slopes, the effective 

reflectivity associated with the edga area la smaller than the reflectivity of 

the center area of the footprint. Therefore, upon reflection from the ocean 

surface, the laser cross-section Is modified by the distribution of surface 
_2 

slopes. The S term in Equation (76) accounts for this modification. The 
mean square width of the received pulse Is obtained from Equation (65) 
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E[o 2 ] " % ♦ \ + ^(tan“ 2 0 T + 2s“ 2 )’ 2 (77) 

c c 

.3 

If 6 Is on ths ordsr of 10 radians or am aller, eha last can In 

1 2 

4* 4 

Equation (75) will ba approximately aqual to — j- tan 0^, Indapandant of 

2 c 

tha saa stats. Therefore, an estimate of 0 ., which Is tha var lanes of 

^ 2 

2 2 4* 4 

surface height can ba obtained, since o^, o f and — tan 6^ ara all 

c 

known parameters of tha altimeter. 

Next, wa taka Into account tha non-Gauss lan nature of the ocean 
surface. Since we are considering normal Incidence, and the divergence 
angle of the laser beam is vary narrow, the specular points that contribute 
to the received field will be those facing upwards or with taro slopes. 
Hence, we can approximate p(£ ,£ ,£) in Equation (23) by 

*• y 

.o - p(o,o,o ( 78 ) 

By using Equation (43) for p(0,0,£), and substituting into Equation (23) for 
Che expected received signal, the integration can be evaluated numerically. 
The results are shown in Figure 8 to Figure 10 for different values of the 
skewness coefficient. 

The mean pulse delay and mean square pulse width in this case are 


_ 2 * . 2 * _ 2 ,2 . 

<T > ■ — + — tan 0_ + — X, o_ 
see T c 3 £ 


E[0 2 ] ■ + a f + ^ 2 ” ®x + ~2 °r^ 1 

c c 


x 2 ) 
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(80) 


where X^ is the skewness coefficient defined previously in Equation (42). 
Jackson [16] pointed out that the value of X^ depends on the wave age 
and Che dominant wave length. Developing seas have a greater skewness 
value, while old seas and swell exhibit smaller-than-cqullibrlum skewness 
values. The equilibrium value of X^ is about 0.2. 
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Figure 9. Mean received waveforms of the laser altimeter for different 
values of the skewness coefficient. 
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Comparing Equation (79) with Equation (76), wa find char* la bias In 
tha aaan of tha altlMtar output. In tanu of altltuda, tha amount of 
blaa la which is about 20X of tha raa vava haight, a significant amount 

consldarlng tha accuracy of tha lasar altlmatar. 

By comparing Equation (SO) with (77) for tha naan squara pulsa width, 

2 

tha dlffaranca is found to ha about 4X of tha Man squara roughnass o f . Wa 
concluda tha actual valua of is irvortant in aatiMtlng tha altituda, 
but is lass critical if only tha saa surface roughnass is to ba axtractad 
f.-on tha racalvad pulsa width. 


37 


6. NON-NORMAL INCIDENCE 

For the analysis in the previous sections, the laser altimeter was 
assured to be pointed at nadir. If the nadir angle is small so that the 
shadowing effects can be neglected, the results of the previous section can 
be easily modified to include the effects of the nadir angle. The system 
geometry is illustrated in Figure 11. The expressions for the means and 
variances of . received pulse delay and width involve 2-D integrations 
over the _? * ( x »y) coordinates, which are transverse to the direction of 
propagation. For the geometry in Figure 11, we have chosen the nadir angle 
so that the propagation axis is normal to the y-axis and intersects the 
x-axis at the angle tt/2 - 4>. In order to apply the results of the previous 
section, we need to determine the apparent altitude z' and apparent surface 
profile £’ for the nadir pointing angle in terms of the actual altitude z 
and profile £. From simple geometrical considerations, we have 
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Again, the divergence angle of the laser beam is small, so that the 
following equation holds for the contributing specular points 


p(S »S) z p(tan $,0.£) 
x y 


( 85 ) 
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Altimeter 



Figure 11. 


Geometry of the laser alcimeter and ocean surface for non-normal 
incidence . 
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By modify log Equation (41) into tha joint danalty function of surfaca 
alopaa and alavatlon avaluatad at ( ■ tan 4, 5 ■ 0, wa hava 

^ y 
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t(2ir 3 ) 1/2 
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( 86 ) 


2 2 S 

In Equation (86) , wa have assuaed that <5^> - <5 y > • — . Using 
Equations (66), (70) and (86), tha expected delay and mean square width 
of the received pulse are calculated 
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( 88 ) 


From Equation (87), we find that the non-zero nadir angle has increased 
the mean delay time by the factor (cos 4) For the mean square pulse 
width in Equation (88) , the last term is the extra spreading due to the tilt 
of the altimeter. The origin of this ten* could be seen from Equation (82), 
where the tilt introduces a linear tern into the effective profile. 

The pulse spreading due to non-noraal incidence will be significant 
if the aagnltude of z tan 0^ tan 4 is comparable to the ras wave height o^, 
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becauae ambiguity will then arlaa In determining from the received pulae 
width. It appear* that the effect of the non-zero nadir angle can be 
minimized by decreaalng the beam divergence angle 6^; however, aa pointed 
out by Gardner [1], thla la achieved at the coat of higher apeckle Induced 


nolae. 
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7. CONCLUSIONS 

In previous chapters, we have shown thet a short pulse laser altimeter 
can be used In the determination of mean sea level and saa stataa. Since 
a much narrower pulse width can be transmitted by a lasar altimeter than by 
a radar altimetar, the Inherent resolution of the laser altimeter is clearly 
higher. We point out that to fully achieve the high accuracy promised by a 
laser altimeter further research on the statistical properties of capillary 
waves and their Interactions with optical radiation has to be done. Also, 
on the true sea surface, whltecaps and foam patches begin to form as wind 
speed increases. Foam patches and whltecaps contain sprays and bubbles 
that cause scattering which Is not easy to analyze theoretically; their 
effects on the received signal may or may not be significant, and should 
be determined by future experimental work. 

Results of Chapter 6 indicate that nadir angle effects enter the 
estimates of both the mean sea level and SWH. When the altimeter is at 
orbital altitudes, it is necessary to measure the nadir angle to a high 
degree of accuracy, since an error on the order of milliradians can cause 
erroneous estimates of the mean sea level and SWH. 
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